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Abstract: 

A microfluidic whispering gallery mode (WGM) biosensing system is pro- 
posed for enhanced delivery and detection of target molecules. A micro- 
toroid resonator coupled to a tapered optical fiber is immersed within a mi- 
crofluidic channel, and supplied with target molecules at various flow rates. 
We show through Monte Carlo simulations that the flow characteristics and 
resonantly enhanced optical forces of the sensor substantially improve both 
the sensing time and yield. When compared to a diffusion-limited sensing 
modality, the average time required to detect a single molecule is reduced 
from more than 100 minutes to less than 10 seconds, and the overall yield 
of the device is enhanced from less than 5% to a maximum of 70.6% for 
femtomolar concentrations of analyte. 
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1. Introduction 

Whispering gallery mode (WGM) resonators are optical sensors that exhibit unprecedented 
sensitivity due to their small mode volume and ultra-high Q factor (Q^ 10 8 ) HHUO. ^ n t ^ le ^ e 
sciences, they have been shown to be excellent candidates for biological sensing 0|3 16), and 
have recently demonstrated single molecule sensitivity by detecting individual proteins Q and 
virus molecules |8| without the use of labels. In a WGM biosensor, light is evanescently coupled 
to a water-immersed resonator via a tapered optical fiber. The resonant frequency of the light 
circumnavigating the resonator is sensitive to small dielectric pertubations in its environment, 
such as the adsorption of proteins on the resonator's surface. Detection of single biomolecules 
is achieved by functionalizing the surface of WGM sensors with anti-bodies or DNA specific 
to the target molecule, and then monitoring the resonance frequency spectrum of the sensor. 

The performance of a biosensor, however, is largely dependent on the transport and binding 
kinetics of target molecules at the sensor's surface. Generally, the transport process at the sur- 
face is exclusively diffusive, and therefore the measurement of sub-picomolar concentrations 
without directed transport can require time scales of hours to days |9|. Introducing a simple 
convective flow over the surface does not appreciably improve the yield of the device, as diffu- 
sion tends to dominate at the sensor surface and convection dominates further away. Since it is 
often desirable to analyze small volumes or extremely dilute concentrations of analyte (9), an 
efficient method of delivering target molecules to the sensor's surface, such as a microfluidic 
system, is needed to ensure an acceptable yield in a practical biosensor. 

In this paper we propose to integrate a microtoroid WGM sensor into a microfluidic system 
that transports molecules directly towards the most sensitve area of the sensor, and then cap- 
tures them by the gradient forces of the resonantly enhanced electromagnetic field. We show 
that these processes combined can overcome the diffusion-limited transport mechanism and re- 
duce the time required to detect femtomolar concentrations of molecules from more than 100 
minutes in a diffusion-limited sensor to less than 10 seconds. As expected, with increasing con- 
vection a compromise is observed between improved sensing time and reduced yield due to 
the shorter time spent by molecules near the sensor surface. However, we show that this com- 
promise between yield and fast detection time can be overcome by designing a microfluidic 
channel geometry that induces vortices for large flow rates, acting to focus the molecules onto 
the sensing region. Significant increases in the yield are then observed over a wide parameter 
regime, with a yield as high as 40% achieved with a 7 second average sensing time. 

2. Theory 

In our proposed sensing system, a molecule of radius a moves toward the WGM sensor with a 
velocity relative to the fluid flow x(t) and experiences a viscous force proportional to its velocity 
by /3 = 6nria, where 77 is the viscosity of the solution. It also experiences a Brownian force 
y/2kbTy(t), where k\ T = 4.1 1 x 10~ 21 J at room temperature, and an attractive optical force 



from the gradient of the evanescently decaying electromagnetic field E(r) within the WGM 
sensor, where a ex is the excess polarizability of the molecule. Thus, the molecule follows a 
trajectory given by 



(1987). 
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where m is the mass of the molecule. When the molecule reaches the sensor, the optical gra- 
dient force F pulls it towards a region of high field intensity. If the molecule adsorbs onto 
the surface, the evanescent field polarizes the molecule, inducing an oscillating dipole moment 
dp = a ex E(r) in the molecule, resulting in a shift of the WGM resonance frequency (Oq given 
by 0: 
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where V m is the volume of the optical mode and E max is the maximum intensity of E(r). The 
minimum frequency shift observable with a given WGM sensor will depend on a range experi- 
mental parameters, including most critically the laser linewidth and power, the WGM Q factor, 
and whether the laser is frequency scanned or locked on resonance. Since detection of sub- 
linewidth resonance shifts is possible when frequency locking and stabilization methods are 
employed II 101 . in the article we select a minimum frequency shift of 1.75 kHz for a detectable 
binding event, roughly 1/10 of the linewidth of a tunable diode laser recently developed by 
Doringshoff et al. ifTTl 

The transport characteristics of the biosensor for a given flow field are determined by numer- 
ically solving Eq. (j2j) to find the expected position of the molecule (x(t)). If the time interval 
dt is small enough that F(x) is approximately constant over the interval and the Brownian term 
(7(f)) =0, a finite difference method provides the next expected position (x n+ \) and velocity 
(x n +i) of the molecule 
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Brownian motion is taken into account by adding Gaussian noise terms with mean zero to the 
position and velocity 1121 . These terms have variances of a 2 = St — j ( \ — exp — ^St^j 

and a, 2 = respectively. When a noise term is added to the position, there is some cor- 
related noise in the velocity given by 8x = 8x+A8x, where 8x is the truly random part of 
the velocity kick and A is the correlation coefficient of the velocity kick in the direction of 
8x. By multiplying the noise in the velocity by 8x and taking the expectation value, we find 
(SxSx) — +A (Sx 2 }. This ordinary differential equation can then be solved to the find 

the correlation coefficient A n+ \ at any step n + 1 given an initial correlation constant A n . 
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Since the variance of the truly random velocity kick (x) is known, we can then add a Brownian 
term to the new velocity (i„+i) given by A„ (8x) +N(0, (x)). 

3. Modeling 

Silica microtoroids were chosen as the WGM resonators that constitute the basis of our biosens- 
ing system, and are fabricated with a combination of photolithography, dry etching and a se- 
lective reflow process [1]. The optical mode volume is much smaller than that of its micro- 
sphere counterpart, and therefore exhibits a higher sensitivity to frequency shifts resulting from 



molecular binding events (Eq. [3]l as well as a larger gradient force (Eq. [TJ. A typical micro- 
toroid is shown in Figure [T]\. We use a finite element model (COMSOL Multiphysics® 3.4) 
to numerically solve for the fundamental transverse magnetic (TM) mode in a water-immersed 
microtoroid with major and minor diameters of D=2Q fim and d=6 flm, respectively. Figure[TJ} 
shows the extent of the evanescent field into the aqueous environment, and illustrates the basis 
for biological sensing using WGM microcavities. To approximate the mode volume of a micro- 
toroid, Spillane et al. [3] derived a simple equation to express the mode volume in terms of its 
"sphere-like" and "step-index fiber-like" components 

-1/3 

(7) 

where Al] 1 '" is the mode area of a step-index fiber with diameter d and Vm Phere is the mode 
volume of sphere with diameter D. Using similiar values as in Ref 13], this expression approx- 
imates the mode volume of a microtoroid with an error less than 4% and, in our case, gives a 
value of V„,=98.4 fiml From Eq. jsj, the maximum frequency shift achievable for a single ad- 
sorbed molecule can then be estimated. We consider the protein Bovine Serum Albumin (BSA), 
with a polarizability a ex = 3.85 x 10~ 27 m 3 . Due to the electromagnetic boundary conditions, 
a discontinuity is present in the electric field at the microtoroid surface, as shown in Figure[T)3, 
enhancing the molecule-resonator interaction. As a result, we find from the numerical solutions 
for the electric field given by COMSOL that E(r) is maximum at the equatorial surface of 
the microtoroid, giving a ratio |E(r)| 2 /|E mav | 2 of one. The maximum frequency shift is then 
predicted to be 6.8 kHz, larger than our minimum detectable shift defined previously. 




Fig. 1. (A) SEM image of a microtoroid. Scale bar is 20 flm. (B) Vertical cross-section of 
the finite element model of the electric field within a a 20 flm diameter toroid immersed 
in water. The circle indicates the boundary of the microtoroid, and the grey area represents 
the supporting silicon pedestals. 
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As shown in Figure |2p, the microfluidic system considered here consists of a 6 flm wide, 2 
flm tall input channel that infuses a 1 1 .6 flm tall chamber which contains the silica microtoroid. 
The input channel has a length of 4 flm. Near the input channel, the chamber is semi-circular 
with a radius of 16 flm, and thus encircles half the microtoroid. At a distance of 16 flm from 
the input channel, the chamber remains 32 flm wide and functions as a large rectangular ouput 
channel with a length of 16 flm. The total volume of the microfluidic system is roughly 6.75 pL. 
COMSOL is used to characterize the fluid flow field within the chamber, with the geometry de- 
fined by a mesh of 58 1 1 points (26305 elements and 3940 boundary elements) and the entrance 
pressure at the input channel defined relative to atmospheric pressure as a boundary setting. 
The flow field around the microtoroid is found by solving the incompressible Navier-Stokes 
equations. Typical horizontal and vertical cross-sections of the flow are shown in Figures |2]\ 
and B, respectively. The flow is quadratic within the input channel, and then spreads out into the 



chamber, giving a low flow field at the surface of the microtoroid and ultimately atmospheric 
pressure at the outlet. 




Fig. 2. (A) Vertical and (B) horizontal cross-sections of the flow field for an input pres- 
sure of 110.2 kPa. White and grey areas indicate the interior of the microtoroid and its 
supporting pedestal, respectively. (C) A microtoroid coupled to a tapered optical fiber and 
incorporated into the microfluidic system. The intensity map around the microtoroid indi- 
cates the intensity of the TM mode at the surface. 

Although the fabrication details of the microfluidic system are not discussed in this paper, 
a flip-chip-based design is proposed as a feasible method to enclose the microtoroid, with the 
microfluidic chamber defined by patterned channels in two identical silicon chips. We envisage 
removing the microtoroid from the top chip via wet etching, leaving the silicon pedestal which 
can then be placed down upon the top of the microtoroid on the bottom chip (Figure |2p). As 
a result, the final device consists of a microtoroid sandwiched by two silicon pedestals (grey 
areas in Figure [2^5). 

To analyze the transport characteristics of our biosensor, we perform Monte Carlo simu- 
lations of BSA proteins passing through the microfluidic system, governed by the equation 
of motion in Eq. Q. A total of 10,000 molecules with initial velocities xq taken to be their 
mean velocities relative to the fluid flow x(t ) are allowed to pass through the system until the 
molecules bind to the microtoroid or flow past it. In the simulation we do not explicity consider 
the density of binding sites (~5 x 10 14 cm~ 2 ) 1 13] on the surface of silica. Rather, we consider 
a molecule bound to the microtoroid when it passes through the mesh boundary defining the 
interface of the chamber and microtoroid, after which it is not allowed to de-bind from the sur- 
face. This is a reasonable approximation in the case of low molecular concentrations, where the 
sensor binding sites are far from saturated. Indeed, our simulations are well within this regime 
when only 10,000 molecules are present. Thus, the kinetics of our binding events are expected 



to be linear, rather than following the typical behavior of first-order Langmuir kinetics (9). A 
time step of 2 /is is used in the finite difference method to track the molecules' positions in time, 
and the frequency shifts due to the binding events are calculated from Eq. ([3]). To investigate 
the effect of convection on binding kinetics we perform a series of simulations in which the 
pressure difference between the input channel and the outlet ranges from to 190 kPa, where a 
difference of 10 kPa roughly corresponds to a flow rate of 0.4 /iL/min. 

For comparison with a diffusion-limited sensing modality, we perform additional simulations 
for a second geometry in which a microtoroid is placed in a micro-droplet of solution with a 
volume Vrf of 1 JJ.L. Since 10,000 molecules suspended in a 1 jxL solution give a concentration 
of about 16 fM, our micro-droplet is representative of a typical miniaturized sensor with respect 
to volume and concentration j9). The volume of the micro-droplet is fixed so that, for compar- 
ison, the time of each diffusive simulation given by V^/q/- is the same as the microfluidic case, 
where q/ is the flow rate of the corresponding microfluidic simulation. Thus, we can compare 
the efficacy of both systems in detecting single molecules for a given sensing time. For both 
geometries, the optical power carried by the tapered optical fiber is varied from mW to 10 
mW, and the operating wavelength chosen for the study is 852.4 nm. 

4. Results 
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Fig. 3. Single molecule binding probability for the microfluidic and diffusion sensing sys- 
tems. The input pressures are (a) 273.2 kPa, (b) 113.2 kPa, (c) 103.2 kPa, (d) 101.4 kPa, 
and (e) corresponds to the response of the diffusion system. The inset shows the transient 
response of the sensor at short time scales. 

In this paper we characterize the efficacy of our WGM bioensor with three properties: the 
yield, or percentage of molecules detected; the efficiency, defined as the probability of detection 
once a molecule reaches the sensor; and the detection time, or the average time needed for a 
molecule to be detected. In Figure [3] we show the number of binding events as a function of 
time, normalized to the total number of molecules for our microfluidic and diffusion sensing 
systems at a fixed laser power of 1 mW. Normalized in this way, the vertical axis corresponds 
to the probability of an initially randomly placed molecule successfully binding to the sensor. 
The average single molecule detection time is found by extrapolating this data to find the time 
at which the binding probability is 50%. 

In the case of diffusion (Figure [3^), a thin depletion zone develops around the surface of the 
sensor and grows as \ft when anisotropy in diffusion is neglected |9|. Initially, the detection at 
the sensor surface also grows as \fi until steady-state accumulation is reached and the detection 



becomes linear with time [14|. By contrast, in the microfluidic system, the sensor response is 
linear with time even for small times and rapidly increases with increasing flow rates. Table [T] 
compares the efficiency and average sensing times of the microfluidic and diffusion systems for 
a laser power of 1 mW. Importantly, the average time required for one molecule to bind with a 
50% probability can be reduced by three orders of magnitude simply by varying the pressure 
difference between the inlet and the outlet of the microfluidic chamber. We note that the asso- 
ciated flow rates are easily attainable in typical experimental conditions using a syringe pump, 
and that the introduction of a modest flow can enable detection of femtomolar concentrations 
on a sub- 100s time scale. Furthermore, over 90% efficiency is maintained for the microfluidic 
system regardless of input pressure, whereas a much lower efficiency (36%) is observed in the 
diffusion system. 

Table 1 . Comparison of the performance of microfludic and diffusion systems as a function 
of input pressure. The incident optical power is 1 mW. 



Pressure difference 
(kPa) 


Efficiency 


Average 
sensing time (s) 


0* 


0.90 


7540 


1.8 


0.92 


671 


11.8 


0.94 


200 


51.8 


0.94 


48 


81.8 


0.94 


18 


171.8 


0.92 


7 


Diffusion Only 


0.36 


6032 



* Atmospheric pressure; transport is exlusively diffusive. 



An important feature of any WGM-based sensor is its ability to attract molecules via the 
gradient of the evanescent optical field. In Table [2j we show that the optical field can signifi- 
cantly reduce the average sensing time of the microfluidic system, and that over 90% efficiency 
is maintained for optical power of 1 mW or more. In Figure |4] we show the effect of increas- 
ing the input laser power on the overall yield of the sensor for a given flow rate and sensing 
time. In the microfluidic system, the presence of 1 mW of incident optical power increases 
the yield from 1 1.2% (0 mW) to 56.7%, illustrating the sheer practicality of employing optical 
cavities as biosensors. It is important to note that most experiments to date HOE I have op- 
erated close to mW, since detection was achieved via repeated scanning of the laser over the 
resonance frequency as opposed to locking to it, resulting in a very low average power within 
the resonator. Increasing the circulating optical power results in a larger evanescent wavefront, 
attracting molecules further away from the sensor throughout the duration of the simulation. 
The yield therefore increases as a function of power, and reaches a maximum of 70.6% for 10 
mW of incident power. By comparison, over the same sensing time scale, less than 5% yield is 
observed when transport is diffusion-limited, regardless of optical power (top of Figure [4]). 

A corollary of the optical force is the redistribution of frequency shifts seen in the histograms 
of Figure [4] A stronger optical field attracts molecules faster, giving them less time to diffuse 
away in the azimuthal direction. Consequently, molecules are more likely to bind to areas of 
the sensor with high field strength. In our simulation, this translates into larger frequency shifts 
from Eq. (BJ, and an increase in the number of counts in the maximum frequency shift bin (6.2 
kHz). The toroidal images in Figure [4] show the final locations of the bound molecules, and 
depict the equatorial localization of bound molecules observed in the simulations. 



Table 2. Performance of the microfluidic system as a function of optical power for a con- 
stant pressure difference of 0.8 kPa. 



Power 


Efficiency 


Average 


(mW) 




sensing time (s) 





0.38 


1421 


1 


0.92 


687 


3 


0.95 


621 


5 


0.96 


607 


10 


0.97 


581 



We find that the microfluidic system always produces a higher yield than its diffusive coun- 
terpart. This can be attributed to the shape of the flow profile after leaving the channel, which 
acts to force molecules onto the sensing area. If the flow rate qf is sufficiently large in the mi- 
crofluidic system, it might be expected that some molecules flow past the sensor faster than they 
are able to diffuse towards it, and therefore fail to bind, leading to a reduction in yield. Figure [5] 
displays the calculated yields for our biosensor as a function of input pressure with mW (yel- 
low squares) and 1 mW (blue circles) of optical power. As expected, the yield sharply declines 
with increasing pressure, since more molecules flow past the microtoroid undetected due to the 
increasing dominance of convection over diffusion. As the pressure is increased above 140 kPa, 
however, vortices develop in the corners of the chamber and input wall, inducing a convective 
flow that directs molecules back towards the front wall of the sensor, as illustrated in Figure [5]3 
and C. Fortuitously the re-directed flow increases the yield of the sensor, which steadily in- 
creases with pressure until a threshold around 275 kPa is reached, after which the flow is again 
less directive and flows past the microtoroid. Thus, careful choice of the microfluidic channel 
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Fig. 4. Comparison of the percentage of molecules detected in the microfluidic and diffu- 
sion systems within 10 seconds. Figures at the top and bottom show the distributions of 
frequency shifts and calculated yields as a function of power for the diffusive and microflu- 
idic systems, respectively. The pressure difference in the microfluidic system is 0.8 kPa, 
and the bin size of the histograms is 2 kHz. 



design allows a very short measurement time to be achieved at high pressure without compro- 
mising the yield of the device. At 273.2 kPa, a detection yield of 40% is achieved with a sensing 
time of 7 seconds (Table [T| for only 1 mW of input power. As can be seen in Figure [5] the op- 
tical force plays less of a role and ceases to increase the yield as fewer stray molecules flow 
past the microtoroid, with the ratios of yields for 1 mW and mW reaching unity at pressures 
greater than ^180 kPa. 

Although the greatest yield of 70.6% is achieved for low convection and a long sensing time, 
the yield associated with high convective flows can be increased by recirculating the flow back 
into the microfiuidic system, so that undetected molecules may be subsequently captured. Thus, 
one would expect an enhancement in the yield increasing with the number of recirculations j, 
allowing for a very high yield Y J to be achieved with a low measurement time: 
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Fig. 5. The effect of pressure on the yield of the microfiuidic system. Insets A-D depict the 
trajectory of the molecules in the microfiuidic chamber, illustrating the rough characteris- 
tics of the flow field within the chamber for the range of flow rates used in this study. The 
incident optical power is mW (yellow squares) or 1 mW (blue circles). 
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where 77 is the efficiency of the sensor, and ^detected and Nj„ are the number of detected and 
incident target molecules, respectively. Note that for an infinite number of recirculations, the 
yield approaches the efficiency as expected, and that Y' is simply the inital yield for one recir- 
culation 0=1). For high convective flows such as the case when the input pressure is 273.2 kPa, 
we calculate that the initial yield of 40% can be increased to more than 75% with only three 
recirculations (and in only 21 seconds). Thus, our microfiuidic system is capable of achieving 
very high yields for very sensible time scales. Furthermore, we note that, even without recir- 
culation, just 1 mW of optical power can increase the yield of the microfiuidic device to more 
than 50% for slow flow rates, more than 20 x greater than the yield from the diffusion-limited 
case. 



5. Conclusion 



We have proposed a WGM biosensing sytem that greatly improves the detection of molecules 
in solution. The chief advantage in using WGM-based resonators as advanced biosensors is that 
single molecules can be detected without the use of labels and with very low input power. Gen- 
erally, the utility of WGM sensors would be limited by the confinement of light to the small vol- 
ume of the resonator. However, we have shown that incorporating a microtoroid WGM sensor 
into a microfluidic system allows for controlled manipulation of molecules and its surround- 
ing fluid, substantially improving the probability of detection over a reasonable time frame. 
Single molecules can be detected with an average measurement time as low as 7 seconds, and 
the overall yield of the device can be enhanced from less than 5% in the case where diffusion 
dominates to as high as 70.6% where the interplay of convection and optical forces becomes 
dominant. Particles or molecules as small as a single BSA protein (roughly 6 nm in radius) can 
be detected, and the analysis can be completed for femtomolar concentrations. These results 
demonstrate that the combination of microfluidics and WGM sensing provides a highly tunable 
system, with both yield for a given concentration and sensing time optimizable by varying the 
input power and flow characteristics of the microfluidic system. 
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